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ABSTRACT 

We present the results of numerical simulations of stationary, spherically outflowing, 
pair winds, with total luminosities of L — 10 35 — 10 42 ergs s -1 . These results have 
direct relevance to the emission from hot, bare, strange stars, which are thought to 
be powerful sources of pairs created by the Coulomb barrier at the quark surface. 
The spectra of emergent photons and pairs are calculated. For L > 2 x 10 35 ergs s _1 , 
photons dominate the emerging emission. As L increases from ~ 10 35 to 10 42 ergs s _1 , 
the mean photon energy decreases from ~ 400 — 500 keV to 40 keV, while the spectrum 
changes in shape from a wide annihilation line to being nearly blackbody with a high 
energy (> 100 keV) tail. Such a correlation of the photon spectrum with the luminosity, 
together with the fact that super-Eddington luminosities can be achieved, might be a 
good observational signature of hot, bare, strange stars. 
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1 INTRODUCTION 

Strange stars made entirely of deconfined quarks have long 
been proposed as a possible alternative to neutron stars 
(e.g., Witten 1984; Alcock, Farhi, & Olinto 1986a; Haensel, 
Zdunik, & Schaeffer 1986). Because strange quark matter 
(SQM) at the surface of a bare strange star is bound by 
strong interactions, not by gravity, such a star can radiate 
with luminosities greatly exceeding the Eddington limit. Lu- 
minosities as high as 10 52 ergs s _1 can be radiated briefly 
after the formation of the star, when the surface tempera- 
ture is Ts ~ 10 11 K (Alcock et al. 1986a; Chmaj, Haensel, 
& Slomihski 1991; Usov 1998, 2001a). Some 10 seconds af- 
ter formation, Ts drops to less than 10 10 K (Page & Usov 
2002). At such temperatures-which may also prevail later 
in the life of the star due to episodes of reheating-electron- 
positron (e ) pairs created by the Coulomb barrier at the 
surface take up most of the thermal emission from the SQM 
surface (Usov 2001a). Some of these pairs annihilate into 
photons, thus forming a pair-photon wind that flows away 
from the star. 

For the typical radius of a strange star (10 6 cm), a 



man 1998; Iwamoto & Takahara 2002). The emerging emis- 
sion consists mostly of photons, so L 7 ~ E. The photon 
spectrum is roughly a black body with a temperature of 
~ 10 10 (E/10 49 ergs s" 1 ) 174 K. The emerging luminosity in 
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tion rate in pairs E ~ 10 43 - 10 49 ergs s" 1 (Usov 2001a) 
For such powerful winds the pair density near the sur- 
face is very high, and the outflowing pairs and photons 
are very nearly in thermal equilibrium almost up to the 
wind photosphere (e.g., Paczyhski 1990). The outflow may 
then be described fairly well by relativistic hydrodynamics 
(Paczyhski 1986, 1990; Goodman 1986; Grimsrud & Wasser- 



L-y ~ 10" 

this applies roug hly down toE~ 10 42 - 10 4 

In contradistinction, for E < 10 42 ergs s" 1 (T s < 9x 10 8 
K), which is the region we explore, the thermalization time 
for the pairs and photons is longer than the escape time, 
and pairs and photons are not in thermal equilibrium. In 
this Letter, we describe (for the first time to our knowledge) 
the results of numerical calculations of the characteristics 
of the emerging emission in pairs and photons in stationary 
winds with energy injection rates E = 10 35 — 10 42 ergs s _1 . 
The details of the set-up, code, processes included, and the 
results for the wind structure, will be described elsewhere 
(Aksenov, Milgrom, & Usov 2003). Here we give a brief ac- 
count. 



2 FORMULATION OF THE PROBLEM 

At the moment of formation of a strange star, the tempera- 
ture in the stellar interior is expected to be a few times 10 11 
K (Haensel, Paczyhski, & Amsterdamski 1991; Cheng & Dai 
2001). The neutrino luminosity of such a hot strange star is 
up to ~ 10 54 ergs s" 1 . The rate of neutrino-induced ejection 
of normal matter, which may be present at the SQM sur- 
face (Alcock et al. 1986a; Glendenning & Weber 1992), is 
very high (e.g., Woosley & Baron 1992; Levinson & Eichler 
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1993). Most probably, in a few seconds after the star forms, 
the normal-matter envelope is blown away completely (Usov 
2001a). The SQM surface remains exposed so as long as the 
surface temperature is higher than ~ 3 x 10 7 K (Usov 1997). 
From the thermal evolution of young bare strange stars it 
follows that this may be the case for at least a few hun- 
dred years after formation, and perhaps a few orders longer, 
depending on the phase of SQM, convection, etc. (Page & 
Usov 2002) 

We consider an pair wind that flows away from a hot, 
bare, unmagnetized, strange star with a radius of R = 10 6 
cm. Pairs are injected from time t — at a constant rate 
into the wind, which is assumed spherical. This results in a 
time-dependent wind which eventually becomes stationary. 

We solve numerically the coupled, relativistic Boltz- 
mann equations for the pairs and photons (e.g., Mezzacappa 
& Bruenn 1993; Aksenov et al. 2003) 

1 dEi fj, d , 2 a , Id 
c at or r dfi 

= £1 

9 



H 2 )f3iEi] 



(1) 



where Ei(e, fi,r,t) = 2ire 3 f3ifi/c 3 are the {r, fi e}-phase- 
space energy densities of e pairs (i = e) and photons 
(i = 7) (fi being the distribution functions); fi = cos 9, 
— 1 < fx ^ 1, where 9 is the angle between the radius- vector 



r and the particle momentum p; /3 7 



1, A= = v e /c 



v/l — (m e c 2 /e e ) 2 ; e is the particle energy; rjt is the emis- 
sion coefficient for the energy production by a particle of 
type i via the physical process labelled by q; and, x\ is the 
corresponding absorption coefficient. The summation runs 
over all considered physical processes. Gravity is neglected 
[but will be included in our more detailed calculations (Ak- 
senov et al. 2003)]. This is a good approximation for the 
energy injection rates we explore, because L ~ 10 35 ergs s _1 
is the Eddington luminosity for the pair plasma, and above 
it, radiation-pressure forces dominate over gravity. 

At the internal boundary, r — R, the input pair number 
flux depends on the temperature Ts at the stellar surface 
alone, and is taken as (Usov 2001a) 
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where < = 20(T S /10 9 K)" 1 . Their energy spectrum is ther- 
mal with temperature Ts, and their angular distribution is 
isotropic for < /1 < 1. The energy injection rate in e ± pairs 
is then 
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where fcs is the Boltzmann constant. 

The stellar surface is assumed to be a perfect mirror 
for both e pairs and photons. At the external boundary 
(r — r cxt ), the pairs and photons escape freely from the 
studied region. 

The thermal emission of pairs used in our simulations 
[eq. (2)] is that expected in the normal, unpaired, phase of 
SQM. However, it has been argued that SQM is a color su- 
perconductor with a critical temperature of ~ 10 11 — 10 12 
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Figure 1. The partial luminosities in pairs (L e ) and photons 
(Lj) as functions of the total luminosity, L. 



K (for a review, see Alford, Bowers, & Rajagopal 2001; Ra- 
jagopal & Wilczek 2001). It is expected that at density much 
higher than the mean nuclear density, p n ~ (2 — 3) x 10 14 
g cm , the color superconductor is in the "Color-Flavor- 
Locked" (CFL) phase, in which quarks of all three flavors 
and three colors are paired in a single condensate. If the 
density is only slightly higher than p n , SQM may be in the 
"2 color-flavor Superconductor" (2SC) phase in which only 
up and down quarks of two colors are paired while the ones 
of the third color and the strange quarks of all three colors 
are unpaired. In the CFL phase, SQM is electrically neutral 
and no electrons are present. So, if at the surface the SQM 
is in the CFL phase, there is no supercritical electric field 
and no pair emission. In the 2SC phase, electrons are 
present, and the thermal emission of e ± pairs from the sur- 
face is practically the same as for normal SQM. Therefore, 
our simulations pertain only to the last two cases. 

Although the injected pair plasma contains no radia- 
tion, as the plasma moves outwards, photons are produced 
by pair annihilation. Other two body processes that we in- 
clude in our simulations are M0ller (e + e + — » e + e + and 
e~e~ — + e~e~) and Bhaba (e + e~ — > e + e~) scattering, 
Comptom scattering (ye — * je), and photon-photon pair 
production (77 — > e + e~). Two body processes do not change 
the total number of particles in a system and thus cannot, 
in themselves, lead to thermal equilibrium. For this rea- 
son, we include bremsstrahlung (ee — * e&y), double Comp- 
ton scattering (■ye — » 767), and three quantum annihila- 
tion (e + e~ — > 777), which change the particle number, even 
though their cross-sections are ~ a -1 ~ 10 2 times smaller 
than those of the two body processes (a is the fine-structure 
constant). 



3 NUMERICAL RESULTS 

We give here the results for the properties of the emerging 
radiation after stationarity is achieved, so the total wind lu- 
minosity is equal to the energy injections rate: L = L e +L 7 = 
E. We present results for different values of L, which is the 
only free parameter. The corresponding surface temperature 
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Figure 2. The energy spectra of emerging photons for differ- 
ent values of the total luminosity, as marked on the curves. Also 
shown, for L = 10 42 ergs s — 1 , is the blackbody spectrum with the 
same energy density as that of the photons at the photosphere of 
the outflowing wind (dashed line). 
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Figure 3. The energy spectra of emerging electrons for different 
values of the total luminosity, as marked on the curves. 



Figure 4. The mean energies of the emerging photons (thick 
solid line) and electrons (thin solid line) as functions of the total 
luminosity. For comparison we show as the dotted line the mean 
energy of blackbody photons for the same energy density as that 
of the photons at the photosphere of the out-flowing wind. Also 
shown as the dashed line is the mean energy of the emerging 
photons in the case when only two body processes are taken into 
account. 

photon energy (e 7 ) from ~ 500 keV at L ~ 10 35 ergs s _1 
to ~ 400 keV at L ~ 10 37 ergs s _1 occurs because of the 
energy transfer from annihilation photons to e pairs via 
Compton scattering (see Fig.0J. As a result of this transfer, 
the emerging e pairs are heated up to the mean energy 
(e e ) ~ 400 keV at L ~ 10 37 ergs s" 1 . For L > 10 37 ergs s" 1 , 
changes in the particle number due to three body processes 
are essential, and their role in thermalization of the outflow- 
ing plasma increases with the increase of L. We see in Figure 
2 that, for L = 10 42 ergs s _1 , the photon spectrum is near 
blackbody, except for the presence of a high-energy tail at 
e 7 > 100 keV. At this luminosity, the mean energy of the 
emerging photons is ~ 40 keV, while the mean energy of the 
blackbody photons is ~ 30 keV (see Fig. 4). (The energy 
resolution in our simulations is ~ 20 keV.) 

For L ~ 10 42 ergs s _1 , the emerging pair energy spec- 
trum is close to a Maxwellian , while for L < 10 40 ergs s _1 
it deviates significantly from it (see Fig. 3) 



Ts is found from equations (2) and (3). Figure shows the 
luminosity fractions in pairs and photons. We see that for 
L > L* ~ 2 x 10 35 ergs s~ , the emerging emission consists 
mostly of photons, while e pairs dominate for L < L*. This 
simply reflects the fact that for L < L* , the pair-annihilation 
time is longer than the escape time, and the injected pairs 
remain mostly intact. At higher luminosities, most pairs an- 
nihilate before escape, and reconversion into pairs is ineffi- 
cient, as the mean energy of photons at the photosphere is 
rather below the pair-creation threshold. 

Figures|5]and|j|]present the energy spectra of the emerg- 
ing photons and pairs for different values of L. At low lumi- 
nosities, L ~ 10 35 — 10 37 ergs s -1 , photons that form in an- 
nihilation of e pairs escape from the vicinity of the strange 
star more or less freely, and the photon spectra resembles 
a very wide annihilation line. The small decrease in mean 



4 DISCUSSION 

We have identified certain characteristics of the expected 
radiation from hot, bare, strange stars that we hope will 
help identify such stars, if they exist. The spectrum, we find, 
is rather hard for the studied luminosity range. This makes 
such stars amenable to detection and study by sensitive, high 
energy instruments, such as INTEGRAL (e.g., Schoenfelder 
2001), which is more sensitive in this range than previous 
detectors. 

As super-Eddington luminosities and, as we find, hard 
X-ray spectra characterize the emission from bare, strange 
stars, soft 7-ray repeaters (SGRs), which are the sources of 
short bursts of hard X-rays with super-Eddington luminosi- 
ties (up to ~ 10 42 -10 45 ergs s 1 ) , are reasonable candidates 
for strange stars (e.g., Alcock, Farhi, & Olinto 1986b; Cheng 
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& Dai 1998; Usov 2001b). The bursting activity of SGRs may 
be explained by fast heating of the stellar surface up to the 
temperature of ~ (1 — 2) x 10 9 K and its subsequent ther- 
mal emission (Usov 2001b,c). The heating mechanism may 
be either impacts of comets onto bare strange stars (Zhang, 
Xu, & Qiao 2000; Usov 2001b) or fast decay of superstrong 
(~ 10 14 - 10 15 G) magnetic fields (Usov 1984; Thompson & 
Duncan 1995; Heyl & Kulkarni 1998). For typical luminosi- 
ties of SGRs (L ~ 10 41 - 10 42 ergs s" 1 ), the mean photons 
energy we find is ~ 40 keV (see Fig. 4), which is consistent 
with observations of SGRs (Hurley 2000). 

Another important idiosyncrasy that we find is a strong 
anti-correlation between spectral hardness and luminosity. 
While at very high luminosities (L > 10 42 — 10 43 ergs s _1 ) 
the spectral temperature increases with luminosity as in 
blackbody radiation, in the range of luminosities we studied, 
where thermal equilibrium is not achieved, the expected cor- 
relation is opposite (see Fig. 4). Such anti-correlations were, 
indeed, observed for SGR 1806-20 and SGR 1900+14 where 
the burst statistic is high enough (e.g., Feroci et al. 2001; 
Gogus et al. 2001; Ibrahim et al. 2001). This is encouraging, 
but a direct comparison of this data with results such as 
ours will require a more detailed analysis. In particular, the 
effects of strong magnetic fields have to be included. The 
observed periodic variations may be due to a rotation of a 
star with non-uniform surface temperature, while our result 
apply for isotropic emission. We hope to deal with this else- 
where. 
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